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ABSTRACT 

Context. The environs of massive, early- type stars have been inspected in recent years in the search for sites where 

particles can be accelerated up to relativist ic energies. Wind regions of massive binaries that collide have already been 

established as sources of high-energy emission; however, there is a different scenario for massive stars where strong 

shocks can also be produced: the bow-shaped region of matter piled up by the action of the stellar strong wind of a 

runaway star interacting with the interstellar medium. 

Aims. We study the bow-shock region produced by a very massive runaway star, BD+43°3654, to look for nonthermal 

radio emission as evidence of a relativistic particle population. 

Methods. We observed the field of BD+43°3654 at two frequencies, 1.42 and 4.86 GHz, with the Very Large Array 

(VLA), and obtained a spectral index map of the radio emission. 

Results. We have detected, for the first time, nonthermal radio emission from the bow shock of a massive runaway star. 

Conclusions. After analyzing the radiative mechanisms that can be at work, we conclude that the region under study 

could produce enough relativistic particles whose radiation might be detectable by forthcoming gamma-ray instruments, 

like CTA North. 
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1. Introduction date close to the O supergiant BD+43°3654 (a,(5[J2000] 



/^ 20^33"^36.077® +43°59'07.4C- / b = 82.41° +2.33°). 

^ ; Early-type stars with high peculiar velocities (i.e. runaway Recently, Comeron & Pasqu'ali (2007) related the star 

^ , stars with velocities ^* > 30 km s , e.g. Gies & Bolton bD+43°3654 to a bow shock detected with the Midcourse 

^ 1986) are uncommon. For instance, Maiz-Apellamz et al. g^^^^ eXperiment (MSX) at D and E bands. They studied 

. (2004, Galactic O star catalog for T/ < 8 stars) lists ^ ^^^ ^^^^^^^ ^^^-^^ ^^^^^-^^ ^^ ^^^ surrounding material and 

^ 8% runaway stars out of 370. These particular stars can proposed the star is a runaway member from Gyg 0B2 

g . be identified by the perturbation they produce m the am- association. Gomeron and Pasquali determined a spectral 

^ ; bient medium (e.g. Kobulnicky et al. 2010 and references ^^^^ 04 If, and derived an age of about 1.6 Myr and a stellar 

. . . therein). When the strong winds of runaway OB stars sweep ^^^^ ^f ^ 70 M©, which makes the star one of the three 

^ ' relatively large amounts of gas and dust, the material piles ^^^^ ^^^^-^^ runaway stars known so far. On the basis 

K^ up in the so-called stellar bow shock. Bow shocks develop ^f ^^^^^ estimates, the authors favor a dynamical ejection 

^^ . as arc-shaped structures, with bows pointing m the same scenario (see Hoogerwerf et al. 2000, 2001 for reviews about 

C^ ■ direction as the stellar velocity, while the star moves su- ^^^ orig;in of runaway stars) 

■ " " ' personically in the surrounding interstellar medium (ISM). ^ j p -d /onno^ ^ • ^ j ^u ^ 

£1, . T n 1 1 .1 r.i To^/r . Gravamazdc & Bomans (2008) suggest instead that 

Ine winds are contined by the ram pressure ot the ISM, at -dt^ , aoo ont:A • ^ r - n \ r -i u i 

,. , p ,, , 1 , -11 . 11 BD+43 3654 is part ot a stellar system, formed by a close 

distances from the star determined by momentum balance. • u • • x- i - • u- • • ^u f 

^, ^ ,, 111 .^1 1- .• 1 . .1 1 . 1 encounter between two tight massive binaries m the core or 

ine stellar and shock-excited radiation heats the dust and ^ r^-oo rr-u ^ u i j u ui ^ i ^ ^ u ^u 

, , ,, , . , ^- - ^ . ^ -. ^ Gyg OB2. 1 he star should be a blue straggler to match the 

gas swept by the bow shock. Ihe dust, m turn, re-radiates ,. i . i j • ^u • u ^u • t? u i • i ^ i 

f, . 1 . r TT-» n timescales involved m their hypothesis. Kobulnicky et al. 

the energy as mid-to-iar IR excess flux. /onin\ j u r ^ • j- i i -^ f rr ^^n a 

^-^ (2010) measured a heliocentric radial velocity or —66. 2 ±9.4 

As soon as IRAS images became available. Van Buren & j^^ g-i^ ^^^ derived a stellar mass-loss rate of 1.6 x 10"^ 

McGray (1988) looked for bow-shaped features near high- y^ ^ vr~^ 

velocity O stars (see al^ Van Buren et al. 1995 and Noriega- \^ analyzed data from the NRAO-VLA Sky Survey 

Crespo et al. 1997). The authors detected an IR candi- (^VSS, Condon et al. 1998). The images revealed a coma- 

shaped source of ~ 7 arcmin, spatially coincident with the 

* Member of GONIGET MSX structure (see Figure 1; NVSS angular resolution: 45" ; 

** Fellow of GONIGET rms noise: 1 mJy beam~^). No point sources above 5<j (40 
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Fig. 1. MSX-D band image (color scale) superposed to 1.4 
GHz-NVSS contours. Levels are: -2, 2 {2a), 5, 8, 11, 15, 
19, 24, 29, 50, 70, and 90 mJy beam-^ 



mJy) that are positionally coincident with the MSX source 
are detected in the MIT-Green Bank Survey (GB6, Griffith 
et al. 1991). Inspection of the continuum emission at 408 
and 1420 MHz with the CGPS Survey (Taylor et al. 2003, 
angular resolutions of 3.4' and 1') confirms that the region 
centered on Cyg 0B2 is complex and has strong emission 
on various angular scales (Peri et al. 2010). 

A radio study of the bow shock can shed light on the 
physical processes that give rise to high-energy emission 
from a stellar source, regardless of the history of the run- 
away star. The shock can accelerate particles up to rela- 
tivistic energies by Fermi mechanism. Energetic electrons 
will cool through synchrotron radiation, producing a non- 
thermal radio source. We carried out radio observations at 
two frequencies to study the nature of the emission from 
the bow shock of BD+43° 3654. 

In this Letter we present the results of the radio observa- 
tions in the form of a spectral index map of the bow-shock 
region of the runaway star and distinguish between ther- 
mal and nonthermal emission regions. We briefly discuss 
the issue of whether a bow shock could produce high-energy 
emission enough to be detected with instruments like Fermi 
or the future Cherenkov Telescope Array (CTA) and the 
conditions that must be fulfilled to achieve detection. 

2. Observations and results 

We carried out continuum observations with the Very Large 
Array (NRAOQ) in two array configurations: C at 1.42 GHz 
in April 2008 and D at 4.86 GHz in August 2008. The set up 
allowed mapping largest structures of the size of the MSX 
source and also guaranteed good matching beams of 12'^ 
The flux calibrator used was 0137+331o. Phase calibrator 
scans -2052+365 at L band and 2007+404 at C band- were 
interleaved with target scans. The total time on source was 
3h at each band and 50 MHz for the total bandwidth. 

The data was calibrated with the AIPS package in the 
standard way and analyzed with the Miriad routines. We 
used imagr to produce robust- weighted images. Figure 2 
presents the resulting images after primary beam correc- 
tion, re-gridded with the same synthesized beam. At both 



^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 

^ Calibrator details can be found at www.vla.nrao.edu/. 
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Fig. 2. Continuum emission at 1.42 GHz (upper panel), 
and at 4.86 GHz (lower panel). Contour levels are —3, 3, 
6, 10, 15, 20, 25, and 60 times the rms of 0.3 and 0.2 mJy 
beam~^. BD+43° 3654 is marked with a cross. The arrow 
represents the velocity of the star, derived from proper mo- 
tions corrected for local motion of the surrounding ISM (see 
text). Synthesized beams of 12'' x 12'' are shown in the top 
right corners. 



frequencies there is emission along the extension of the 
MSX source. The radio source is larger at 1.42 GHz, toward 
the eastern half of IR contours of the bow shock (increasing 
right ascensions). Clearly, the spectral index a {S^ oc u^) 
is negative in there. In both images, the rms attained is 
similar (0.3 mJy beam"^ at 1.42 GHz and 0.2 at 4.8, 12" 
synthesized beam). 

A detached ellipsoidal source at Qf,J[J2000] = 20^ 
33"'35^44°04'30", called ES here, is the strongest one in 
the field. The fluxes are 6'es(1.4GHz) = 105 ± 10 mJy and 
6'es(4.8GHz) = 95 + 5 mJy. ES spectral index < a >= 
—0.1 + 0.1 is characteristic of an optically thin Hll region. 

The hypothesis of a physical association between the 
star and the radio/IR features is supported by the very 
good agreement of the residual proper motion of the star 
(15.7° + 9.4° east of north, see Comeron & Pasquah 2007) 
and the direction from the star to the apsis of the bow shock 
(8.8° + 10° east of north). The star velocity vector on the 
plane of the sky is represented in the upper panel of Fig. 2. 

We used the continuum images at 1.42 and 4.86 GHz to 
build a spectral index distribution map. We only considered 



Benaglia et al.: Nonthermal emission from the bow shock of a runaway star 







1 , , , 


, , , 


II- 


4°05' 




\ 




- 


+ 04' 








- 


+ 03' 


. 


i^^^i 


^\ 


- 




4 


^ 


4^^ 


\ 


+ 02' 


m 


"~~^ 




I 


+ 01' 


^ , 1 , , 




, , , 


- 



- J-0.20 H 



20^^33^50' 40' 30' 20' 

Right Ascension (J2000) 



\ 



$ 




} 



0.060 >< 



\ 



Right Ascension (J2000) 



Fig. 3. Top: Spectral index distribution. Bottom: Spectral 
index error distribution. 
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Fig. 4. Acceleration ('Ace'), escape ('Esc'), and cool- 
ing times for electrons, due to synchrotron radiation 
('Synchr'), due to inverse Compton scattering of dust pho- 
tons ('IC (dust)'), stellar photons ('IC (star)'), and back- 
groud photons ('IC (bg)'). Coohng time for relativistic 
Bremsstrahlung radiation indicated as 'Rel B'. The figure 
is for the case with equal energy density in electrons and 
protons (a = 1, see text). 



input pixels with a signal-to-noise ratio > 4. Besides this, 
the spectral index map was masked for a signal-to-noise 
ratio > 10. Figure 3 shows the spectral index distribution 
and corresponding noise maps. 
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Fig. 5. Spectral energy distribution for the case where the 
density of primary electrons is equal to the density of pro- 
tons (a = 1). Acronyms as in Figure 4. Measured radio 
fluxes from VLA observations: 'VLA Obs'. The MSX lumi- 
nosity at D band is also represented. The contribution from 
secondary pairs is negligible in this case, so is not shown 
here. 



3. Nonthermal radio emission from the bow shock 

In Figure 3 we have represented the spectral index values 
derived from 1.42 and 4.86 GHz data at positions where 
(i) the radio continuum signals are well above the noise (4 
times or more) and (ii) the values of the spectral index er- 
rors are low (< 10%). It can be appreciated from the figure 
that there is a significant fraction of the portrayed area that 
shows negative spectral indices. The a- values range from 
—0.8 to +0.3, approximately from W to E. The source ES, 
in contrast, is characterized by a ^ 0, and it looks unre- 
lated to the bow-shock feature. 

In what follows, we assume that the MSX-D&E bands 
excess emission represents the projection of the bow shock 
in the plane of the sky. We adopt a distance to the bow 
shock of 1.4 kpc (Hanson 2003, Comeron & Pasquali 2007). 
We measured a distance from the star BD+43°3654 to the 
bow shock i? = 5^ or 2 pc (at 1.4 kpc). The shape of 
the bow shock can be fitted to a spherical cap of height 
2' = 0.8 pc. We approximated the area of the MSX struc- 
ture on the plane of the sky by 6' in length times 2' in 
height. If the depth is equal to the length, the volume oc- 
cupied by the bow shock is ~ 4.6 pc^. Kobulnicky et al. 
(2010) suggest a conservative lower limit for the particle 
density in the bow-shock region of 100 cm~^. The area cir- 
cumscribed by the projection of the cap is ~1.8 pc^. Most 
of this area shows a source of nonthermal radiation with 
averaged index a = —0.5. Such a radiation is expected 
from synchrotron emission generated by relativistic elec- 
trons accelerated either at the forward shock in the ISM 
or in the reverse shock in the stellar wind. Since the lat- 
ter is expected to have velocities ^rs ~ ^w ^ 2300 km s~^, 
it should be more efficient for electron acceleration (as the 
efficiency rj oc (vrs/c)'^ for first-order Fermi-like diffusive ac- 
celeration, e.g. Drury 1983). The forward shock moves in a 
denser medium with lower velocities, ^ 400 km s~^, but it 
might be effective to accelerate electrons as well although 
only up to lower energies. 
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To estimate the number density of relativist ic particles, 
we followed the standard method by using the observed 
flux density and spectral slope, along with the hypothe- 
sis of equipartition between magnetic fleld and relativis- 
tic particles energy density (e.g. Guinzburg & Syrovatskii 
1964, Araudo et al. 2007). We considered that the energy 
density of relativistic particles makes three contributions 
(addends): 



Eeinei{Eei)dE, 



el 



- / Epnp{Ep)dEp + / Ee2ne2{Ee2)dEe2, 



(1) 



where el, p, and e2 stand for relativistic primary electrons, 
protons, and secondary electron-positron pairs (i.e. pairs 
coming from charged pion decays), respectively, and n is 
the number density. The relation between primary electrons 
and protons is Up = auei^ with a > 0. Three cases were con- 
sidered: a = (just electrons), a = 1 (equal energy density 
in both species), and a = 100 (proton-dominated case, as 
observed in the galactic cosmic rays). The order of magni- 
tude of the equipartition magnetic fleld led to5~5xlO~^ 
G (three times this value for the case a = 100). The max- 
imum value for the energy of the particles was determined 
through the balance of energy gain and losses. Different loss 
mechanisms were considered: synchrotron radiation, inverse 
Compton (IC) scattering of IR, stellar, cosmic microwave 
background photons, relativistic Bremsstrahlung, and par- 
ticle escape from the radiation region due to convection by 
the stellar wind. In the case of protons, the only relevant 
losses are proton-proton {pp) inelastic collisions and con- 
vective escape. Diffusion is negligible in comparison to con- 
vection in this situation. Both primary electrons and pro- 
tons reach energies up to ^ 10^^ eV, which is imposed by 
nonradiative losses, except for a = 100, where synchrotron 
losses dominate for electrons. In Fig. |3]we show the losses 
for electrons in the case a = 1 . Values of magnetic fleld and 
maximum energies obtained for electrons and protons are 
given in Table 1. 



Table 1. Magnetic fleld and maximum energies calculated 
for the three cases of a (see text) for primary electrons and 
maximum energy for protons. 



B 



:V] 



"O 4.0 X IQ-^ 8 X 10'^ 



eV] 



1 5.0 X 10"^ 
100 1.4 X 10"^ 



9 X 10^^ 9 X 10^^ 
2 X 10^^ 3 X 10^^ 



a. Detailed calculations for a set of main parameters will 
be given elsewhere. The pp contribution extends well into 
the TeV regime, but it is weaker and will be difflcult to 
detect with ground-based Cherenkov telescope arrays like 
VERITAS or MAGIC II. In contrast, if the relativistic par- 
ticle contain is proton-dominated (a = 100), gamma-rays 
from neutral pion decay dominate the high-energy spec- 
trum. The planned CTA North observatory might detect 
the source, easily yielding imformation on the cutoff at high 
energies. Observations of the spectral slope in this regime 
can be used to identify the proton content through the lu- 
minosity level and the proton spectral index, since it is 
preserved in the corresponding photon index. Radio po- 
larization data will provide additional information on the 
magnetic fleld. Observations of BD+43° 3654 with X-ray 
observatories like XMM-Newton and Chandra can be very 
useful for determining the cutoff of the synchrotron spec- 
trum, which is directly related to the maximum energy of 
the electrons. This, in turn, would yield valuable informa- 
tion on the actual value of the magnetic fleld and the cor- 
rectness of the equipartition hypothesis. 
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4. Discussion and perspectives 

The presence of highly relativistic particles in a dense 
medium with high photon density can result in the effl- 
cient generation of gamma-rays. Although protons can be 
effectively accelerated up to the highest energies only in the 
shocked wind, where the density is low, they can diffuse or 
be convected upstream up to the region with the swept ma- 
terial and densities of n ~ 100 cm~^. The corresponding 
gamma-ray emissivity can be calculated using the delta- 
functional approximation (e.g. Aharonian & Atoyan 2000, 
Kelner et al. 2006). For the case a = 1, the total luminosity 
from pp interactions is similar to what is obtained from rela- 
tivistic Bremsstrahlung of electrons, since the cross sections 
are similar. In Figure [5] we show the spectral energy distri- 
bution obtained for the case a = 1, with all contributions 
included (synchrotron self-Compton losses are negligible). 
It can be seen that in this case the inverse Compton up- 
scattering of IR photons is the major contribution at high 
energies, with a peak around 100 GeV. The detectability 
of the source by instruments like the Fermi 7-ray observa- 
tory LAT will depend on the actual particle density and 
the contribution related to the secondary electrons at large 
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